(19) 



J 



Europaisches Patentamt 
European Patent Office 
Office europeen des brevets 




(12) 



(43) Date of publication: 

14.06.2000 Bulletin 2000/24 



(n) EP 1 009 124 A2 

EUROPEAN PATENT APPLICATION 

(51) lntCl7: H04L 1/06 



(21) Application number: 99309613.0 

(22) Date of filing: 30.11.1999 



(84) Designated Contracting States: 

AT BE CH CY DE DK ES Fl FR GB GR IE IT LI LU 
MC NL PT SE 

Designated Extension States: 
AL LT LV MK RO SI 

(30) Priority: 07.12.1998 US 206843 

(71) Applicant: LUCENT TECHNOLOGIES INC. 
Murray Hill, New Jersey 07974-0636 (US) 

(72) Inventors: 

• Hochwald, Bertrand M. 
Summit, New Jersey 07901 (US) 



• Marzetta, Thomas Louis 
Summit, New Jersey 07901 (US) 

• Richardson, Thomas J. 

South Orange, New Jersey 07079 (US) 

• Sweldens, Wim 

New Providence, New Jersey 07974 (US) 

• Urbanke, Rudiger L. 

Murray Hill, New Jersey 07974 (US) 

(74) Representative: 

Buckley, Christopher Simon Thirsk et al 

Lucent Technologies (UK) Ltd, 

5 Mornington Road 

Woodford Green, Essex IG8 0TU (GB) 



(54) Wireless transmission method for antenna arrays ftsjng^ji^iy^pace^time signals 



(57) Messages according to a ■ynethodiof^wineless 
signaPtransmission are sent as sequences -W-uT^itary' 
spacenime^signals. Eacbisjgnal is representable as a 
matrix in which each column represents a respective 
transmitting antenna of a transmitting array, and each 
row represents a time segment. Each such matrix is pro- 
portional to a complex-valued matrix, all columns¥of„B 



v^Jcl&ai&JTHfii^^ A new method for cre- 

ating a cojistjBjJ^ti^ comprises pro- 

viding an initial signal in the form of a complex, orthon- 
ormal matrix. The method further comprises generating 
a further plurality of matrices by a process that assures 
that each of the generated matrices is related to the in- 
itial matrix as a product of one or more multiplications 
of the initial matrix by one or more unitary matrices. 
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Description 

Fieid of the Invention 

5 [0001] This invention relates to modulation methods for wireless signal transmission. More particularly, the invention 
relates to modulation methods that reduce the error rates of received signals in fading environments and that enable 
data rates to be increased without the need to increase bandwidth or transmitted power. Still more particularly, the 
invention relates to such methods in conjunction with the use of multiple antenna arrays. 

10 Background of the Invention 

[0002] It is generally desirable to reduce error rates, and to increase transmission rates, in wireless transmission 
systems. Multiple-antenna arrays can be used to achieve these desirable effects. 

[0003] Fading is one of several physical phenomena that tend to increase error rates, or to reduce channel capacity, 
15 in wireless transmission systems. Fading is the result of destructive interference, alt the receiver, between correlated 
signal portions that because of scattering have arrived over different-length paths. 

[0004] One technique that tends to mitigate the effects of fading is differential phase modulation, in which phase 
differences carry transmitted information. Although differential phase modulation is a known technique for single-an- 
tenna transmission and reception in fading environments, there are no known adaptations of this technique for use 

20 with multiple-antenna arrays. 

[0005] However, in certain fading environments, the theoretical capacity of a multiple-antenna communication link 
increases linearly with the size of the transmitter or receiver array, this effect being determined by the array having the 
lesser number of antennas. This effect has been predicted for rich scattering environments in which fading is "flat." 
That is, the propagation coefficients that describe the effect of the physical transmission channel on the transmitted 

25 signal are approximately independent of frequency over the signal bandwidth. Flat fading can be achieved in practice 
for a particular environment if the bandwidth is not too great, or if it is restricted appropriately. 

[0006] Significantly, such a linear increase in capacity occurs only if the propagation coefficients between all pairs 
of transmitter and receiver antennas are known to the receiver. In practice, this condition can be met only if the receiver 
is trained, from time to time, by receiving known training signals from the transmitter. 

30 [0007] Communication methods that use such a training procedure are described, for example, in the co-pending 
U.S. Patent Application Serial No. 08/938,168, commonly assigned herewith, filed on September 26, 1997 by B.M. 
Hochwald et al. under the title, "Multiple Antenna Communication System and Method Thereof." 
[0008] Other co-pending patent applications, commonly assigned herewith, that describe related subject matter are 
Serial No. 08/673,981 , filed on July 1 , 1 996 by G.J. Foschini under the title "Wireless Communications System Having 

35 a Layered Space-Time Architecture Employing Multi-Element Antennas, " Serial No. 09/060,657, filed on April 1 5, 1 998 
by G.J. Foschini and G.D. Golden under the title "Wireless Communications System Having a Space-Time Architecture 
Employing Multi-Element Antennas at Both the Transmitter and Receiver," and a patent application filed on July 10, 
1998 by T.L Marzetta under the title "Determining Channel Characteristics in a Space-Time Architecture Wireless 
Communication System Having Multi-Element Antennas." 

40 [0009] Unfortunately, training intervals cut into the available time during which data may be transmitted. The length 
of this interval increases as the number of transmitter antennas is increased. Moreover, the propagation coefficients 
can be treated as constant only over an average period of time referred to as the fading coherence interval. To be 
effective, training should be repeated at least once per such interval. However, fading is very rapid in some environ- 
ments, such as those in which a mobile station is operating within a rapidly moving vehicle. For rapid fading environ- 

45 ments, the time between fades may be too short for the communication system to learn the propagation coefficients 
belonging to even one transmitting antenna, much less those of a multiple antenna array 

[001 0] Thus, there remains a need to more fully realize, in practice, the theoretical benefits of multiple antenna arrays 
in fading environments. 

[0011] In the co-pending U.S. Patent Application Serial No. 09/134,297, commonly assigned herewith, filed on August 
50 1 4, 1 998 by B.M. Hochwald et al. under the title, "Wireless Transmission Method for Antenna Arrays, Having Improved 
Resistance to Fading," there was described a new method of signal modulation. This new method, which we refer to 
as "Unitary Space-Time Modulation (USTM)," is robust against fading and receiver-induced noise in flat fading envi- 
ronments. Significantly, it does not require knowledge of the propagation coefficients, although in some implementa- 
tions, such knowledge can be used to further improve performance. 
55 [001 2] In USTM, each message to be transmitted is transformed into a sequence of signals selected from a constel- 
lation of L possible signals, L a positive integer. (Thus, each transmitted signal embodies a number of bits given by 
log L. In the present discussion, "log" will denote the binary logarithm.) Each of these symbols is, itself, a time sequence 
of complex amplitudes for transmission by the transmitting antenna or antennas. (We will speak, in general terms, of 
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-atransm.tung array having a plural transmitting antennas. However, it should be noted that the number M of transmitting 
antennas may be 1.) The transmissions by all of the antennas in the transmitting array are concerted All of these 
transmissions (for a given signal) are made in'the^same sequence of T successive*™ units (which we refer to as 
symbol intervals), T a positive integer. ^ * 

5 [001 3] Thus, a signal may be represented by a complex-valued matrix having T rows and M columns Each column 
corresponds to a respective antenna of the transmitting array, and represents the sequence of complex amplitudes to 
be transmitted by that antenna. Each row corresponds to a particular one of the T symbol intervals, and describes the 
complex amplitude to be transmitted by each respective antenna during that interval. Such a set of complex amplitudes 
is reterredto as a "symbol.- Each.symboAis distributecnnspace (i.e.. across the transmitting array), and^achs/ona/ 

'o is ccimposed:of:T£symboJs:dist^buted in time. y 
[0014] Significantly, each signal matrix must have the property that all of its columns are orthonormal (It should be 
noted ,n th.s regard that corresponding to a signal matrix <D. the baseband signals provided to the transmitting array 
are represented by matrix S. where S = 7tP4>. Here. P is the average power fed into each antenna.) Because each 

« ™< c? 6 C T °J UmnS 9,h X th6re ° an n8Ver be more ,han T such columns tnat are a" orthogonal to each other 

[001 5] There are L signals, and M columns per signal. Thus, over the entire constellation, there are LxM columns 
Because there will typically be many signals in the constellation (constellation sizes in the hundreds of thousands or 
even more, are des,rable in at least some applications), LXM will typically be much greater than T. Well known math- 
ematical properties dictate that there can be no more than T mutually orthonormal column vectors. Therefore it will 
be unlikely that, given a randomly chosen pair of signal matrices, the columns of one such matrix will be orthogonal to 

* u the columns of the other. 

[001 6] If such orthogonality between the respective columns of signal pairs were possible, the probability of confusing 
one received signal for another would be reduced to its ideal minimum value. Given that this ideal condition is unat- 
tainable, it is desirable, instead, to design the signal constellation in such a way that correlations between pairs of 
signal matrices, of a kind that tends to increase the error probability, are made as small as possible 
[0017] U.S. Patent Application Serial No. 09/1 34,297, cited above, describes techniques for minimizing these cor- 
relations that are most useful when the number M of transmitting antennas is relatively small. What has been lacking 
until now, is a more powerful technique that can readily generate signal constellations of low correlation when M L 
and T are relatively large, without demanding impractical amounts of computational resources. 

30 Summary of the Invention 

[0018] Such a solution is described here. 

[0019] When a signal matrix (of dimensionality TXM) is left-multiplied by a TXT unitary matrix e, the product is a 
new TXM matrix, which also has the column-orthonormality properties that qualify it to serve as a signal matrix In 
general, the new TXM matrix can again be left-multiplied by ©, and so on ad infinitem, to generate many TXM matrices 
having orthonormal columns. 
[0020] *fe*xm*&mwetoh^ 
«9^f^tivelyitow-correla^ 

s^tl&gsii^yf ^ — — 

40 f? 0211 / particular class of TXT unitar V matrices that we have found useful in this regard is the class of matrices 
that are/th roots of the TXT identity matrix. It will be appreciated that if 0 belongs to such a class, no more thani 
distinct signal matrices can be generated from an initial signal matrix by repeated applications of 0. The parameter I 
may be chosen to equal the size L of the desired signal constellation, or it may be chosen to be larger than L If i is 
larger than L, the constellation is assembled by selecting a subset of size L from the generated matrices 
rj0022] One useful calculation of a pairwise correlation between signal matrices <D /a nd <*,, involves taking the product 
mn£' S *Sf rS *!] e su P erscri P l f deno,es «he conjugate transpose. We refer to this product as the correlation matrix. 
[0023] The indices /and /'are used to denote particular signal matrices. Thus, each of these indices takes on integer 
values from 1 to L. Suppose now that the initial signal matrix is indexed by 1. the result of one multiplication by 0 is 
indexed by 2, and so on. Thus, the signal matrix <b, is the result of M multiplications by 0. Then if © belongs to the 
class of ^th roots of the identity matrix, the correlation matrices will have the following property: Given indices a B y 8 
[a-B]mod^ = [ Y - 5] mod J, then «I>t o 0(J = 0 t <t> s . When this property holds, we say that the signal matrices'have a 
circulant correlation structure. 
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[0024] Accordingly, the invention in one aspect involves a^rjejhodjojicr.eating a constellation!^ 
f Ua ^ m,ssron T*!^^ 

|jwhich are orthonormal. The method further comprises generating a further plurality of matrices by a process that as- 
sures that each of the generated matrices is related to the initial matrix as a product of one or more multiplications of 
|the initial matrix by a unitary matrix. 

10025] In another aspect, the invention involves a method for the wireless transmission of signals. According to such 
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method, at least one baseband signal is generated, and then it is transmitted, on a radio-frequency carrier, from an 
array of one or more transmission antennas. Each transmission of a baseband signal is carried out by transmitting a 
sequence of amplitudes from each of the one or more antennas of the array. 

[0026] Each such transmission is made in accordance with a signal matrix having columns and rows : in which each 
5 column represents a distinct antenna of the array, and each row represents a distinct time interval, such that the entries 
in each column are proportional to the amplitudes to be transmitted, in sequence, from a corresponding antenna of 
the array. The columns of each signal matrix are orthonormal. 

[0027] Each signal matrix is selected from a constellation of available signal matrices. Each matrix in the constellation 
is either an initial matrix, or it is related to the initial matrix as a product of one or more multiplications of the initial matrix 
io by a unitary matrix. 

[0028] In alternate embodiments of the invention, the matrices of the signal constellation are related to the initial 
matrix as products of multiplications of the initial matrix by one or more unitary matrices. 

[0029] In particular embodiments of the invention, the signal constellation is a subset of a set of matrices having a 
circulant correlation structure. 

75 

Brief Description of the Drawing 

[0030] FIG. 1 is a schematic block diagram of a communication system for transmitting and receiving signals in 
accordance with the invention. 

20 [0031] FIG. 2 is a flowchart of a method for generating a signal constellation, according to the invention in an illus- 
trative embodiment. 

[0032] FIG. 3 is a high-level flowchart of an alternate method for generating a signal constellation. 

[0033] FIGS. 4 and 5 are graphs of the results of simulation experiments, illustrating predicted error rates using an 

exemplary signal constellation created in accordance with the invention. 

25 

Detailed Description 

[0034] FIG. 1 shows a baseband signal 10 input to a transmitting array of antennas 15.1-15.3, and transmitted to a 
receiving array of antennas 20.1 , 20.2. Thus, in the communication system shown, M = 3 and N = 2. (N is the number 

30 of antennas in the receiving array, which may be 1 or greater.) It should be noted that although one array is identified 
here as transmitting and the other as receiving, the principles of the invention will apply to bidirectional as well as to 
unidirectional communication systems. The physical transmission channel between the transmitting and receiving an- 
tennas is characterized by a set of M-N propagation coefficients h^, i = 1,..., M, y = 1,..., N, each a complex scalar 
characterizing the response at receiving antenna /due to transmissions from transmitting antenna /. 

35 [0035] At each value r of discrete time, r = 1 T, one of the rows of the signal matrix is input to the transmitting 
array. FIG. 1 shows the fth such row being input, with each entry in the row input to a respective one of antennas 
1 5. 1 -1 5. 3. Each entry of the signal matrix represents a complex-valued baseband voltage level which, for transmission, 
is modulated onto the carrier frequency according to known methods. 

[0036] At each receiving antenna 20. 1 , 20.2, the antenna response is amplified and demodulated to baseband ac- 
40 cording to known methods. Receiver noise, which is here assumed to be statistically independent among the N receivers 
and T symbol periods, is represented in the figure as a component w n added to the output of antenna 20.1 and a 
component w a added to the output of antenna 20.2 at each time r. After demodulation to baseband, the output of the 
antenna array at time t is where n = 1 for antenna 20.1 and n- 2 for antenna 20.2. In vector notation, the response 
X, of the receiving array to the fth row S f of transmitted signal matrix S (the index / is suppressed here) is given by X, 
45 = S f H + If H can be treated as constant during the time period T, then over that period, the response of the receiver 
array is given by X = S H + W, where W is a TXN matrix whose t,n entry represents the additive noise at time f and 
receiver n. 

[0037] An important consideration in any method of message transmission is capacity; that is, the amount of infor- 
mation that can be reliably transmitted per unit time in a communication channel. We define a channel use as one 
so block of T transmitted symbols (i.e., one transmitted signal matrix). We measure the channel capacity in bits per channel 
use. Data can be transmitted reliably at any rate less than the channel capacity. Thus, the channel capacity limits the 
number of bits per signal, or log L. 

[0038] As noted, the coherence intervals the length of time over which the propagation coefficients can be regarded 
as approximately constant. In the following discussion, the symbol x will represent coherence interval. 
55 [0039] The number T of symbols per signal should not exceed t, because if it does, fading effects will tend to corrupt 
the received signal, and the advantages achievable with our method will be reduced. However, if T is substantially less 
than x, the channel will be used inefficiently, which might lead to a greater error rate. Therefore, it will generally be 
advantageous for T to be equal, or nearly equal, to x. 



SDCOD: <EP 10091 24A2 I > 



4 



10 



IS 



20 



25 



30 



3S 



40 



45 



SO 



SS 



EP 1 009 124 A2 

[0040] We performed a theoretical analysis of channel capacity achievable with our method, in a fading environment 

noV^oo^Tl ° f a i!!? Ve . Ce,Ver n ° iSe - ThiS analySiS is discussed in ^e co-pending application Serial No. 
09/134,297, cited above. Although this analysis was based upon a theoretical model that contained a number of sim- 
plifying assumptions, it led to two predictions that we believe will be valid for many transmission systems of practical 
interest. These predictions are. 

1) Incasing the^umbemof lransmW 

2) There are two .ndependent conditions under which the transmission rate obtainable using USTM will be no 
worse than that obtainable using any other modulation method: 

(a) When T » M . 

(b) When T > M and the signal-to-noise ratio is sufficiently high. 

[0041 ] The received signal X may differ in dimensionality from the transmitted signal. Moreover, the received signal 
will generally show the effects of uncertainty contributed by the propagation coefficients and the additive receiver noise 

, '! advanta 9 eous to use a decj s'on Process at the receiving end to restore the transmitted signal 
[0042] The decision processes known in the art as "maximum likelihood (ML) receivers," are useful in this regard A 
ML receiver functions by selecting that candidate signal that maximizes the likelihood of observing the signal X actually 
received. A second decision process, also known in the art. is a maximum a posterior/probability (MAP) receiver If all 
of the candidate signals are equally likely, the ML and MAP receivers are equivalent. If the candidate signals are not 
equally likely, it is advantageous to employ the MAP receiver, which is readily achieved through a simple and well- 
known modification to the ML receiver. 

[0043] With either receiver, a score is computed for each candidate signal, in accordance with a decision process 
The received signal is then identified with candidate signal that yields the greatest score 

that.was transmrtted,Th,s procedure comprises calculating the conditional probability p(X !<[>,) of receiving this particular 
X, given that the transmitted signal matrix was each of the <*, in turn. The 0,that yields the greatest value of this 
conditional probability is identified as the transmitted signal. This "maximum likelihood signal" o^- js symbolically 
represented as the argument of a maximization procedure by the expression 

tfc* 1 - = argmax p(K I *,). 

. ^i!!^ f ^P^ation coefficients are Rayleigh distributed, theimaximum likelihood signal can be evaluated : 
by.max.mizmg ^ arpamcuMrly^imple^^sibn: (By "Rayleigh distributed," we mean that: the values oflhe propagation 
coefficients n mn are identically distributed and statistically independent, that the magnitudes of the propagation coef- 
ficients are proportional to Rayleigh distributed random variables, and that their phases are uniformly distributed from 
0 to i 2k radians. Those skilled in the art will understand from this that the real and imaginary parts of the propagation 
coefficients are zero-mean, independent, identically distributed, Gaussian random variables ) 
[0046] In such a case, 

M n 



argmax p(X | 4>,) = argmax £ J)|*LiC 



In the double-sum expression, the quantity within the vertical bars is ihe vector dot product between the complex 
conjugate of the nrth column of <J>, and the nth column of X. (The underscore beneath the factors within the vertical 
oars denotes that these quantities are column vectors.) 

[0047] in our theoretical model, we have been able to find a convenient upper bound, referred to as a Chernoff upper 
STt-S *S ,wo " S '5T 1 error Probability; i.e., for the probability of error given two signals <t>„ 4>, transmitted with equal 
probability. We found that this bound depends only on M, T. N, the signal-to-noise ratio p. and on the M quantities d 
that are the singular values of the MxM correlation matrix <b]t <t>,. 

[0048] The "singular values" are understood as follows. It is a basic result from linear algebra that any matrix can 
be factored into the product ABCt , j n which A and C are unitary matrices, and B (not necessarily a square matrix) 
has the following properties: all non-diagonal entries are 0, all diagonal entries are real and non-negative, and the 
diagonal entries occur in decreasing order. The diagonal entries are the singular values of the original matrix 
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[0049] Broadly speaking, the singular values of <J>,t 0> r are a measure of the similarity of the sub-spaces spanned by 
the columns of the respective matrices. That is, linear combinations of the M columns of the unprimed matrix generate 
an M-dimensional subspace of T-dimensional complex space. Similarly, linear combinations of the M columns of the 
primed matrix generate a different M-dimensional subspace of T-dimensional complex space. For distinct transmitted 
signals to be discriminated readily and with high assurance, these respective subspaces should be, in a certain math- 
ematical sense, as dissimilar as possible. The smaller the singular values, the greater this dissimilarity. Thus, decreas- 
ing any given singular value (with, e.g. , the others held constant) will tend to decrease the probability of error in decoding 
the received signal. Methods for obtaining the singular values of a matrix are well known, and need not be described 
here. 

[0050] The Chemoff upper bound CU B. is expressed by 



where p represents the signal-to-noise ratio (in units of power/power), and the other symbols are as defined above. 
According to our theoretical model, the two-signal error probability will never be greater than this quantity For a given 
signal-to-noise ratio and a given set of choices for T and M, this quantity depends upon the singular values d m 
[0051] The C.U.B., and thus the maximum possible error probability, is minimized when all of the singular values are 
driven as low as possible, preferably to zero. However, it will not generally be possible to drive more than a few, if 
indeed any, of the singular values to zero, over all pairs of matrices <X>,, O,. . That is, all of the singular values will be 
zero only if all M columns of each signal matrix of the signal constellation are orthogonal to all M columns of every 
other matrix in the constellation. This, as noted, is not generally possible. Therefore, the design of signal constellations 
will generally be aimed at producing a set of singular values that, according to an appropriate measure, are collectively 
as small as possible. 

[0052] We have devised a figure of merit, which we denote by the symbol 8, which expresses, in a collective sense, 
how small the singular values are, over all signal pairs taken from a constellation. For a given correlation matrix O,* 
Q> h define the norm 



Then 5 is the greatest value taken by this norm, over all possible pairs of signals in the constellation. That is, 



[0053] It should be noted that the probability of error is invariant to transformations of the signal constellation that 
are carried out by left-multiplying all signals in the constellation by a common, TXT unitary matrix. (Left-multiplication 
is sometimes referred to as "pre-multiplication.") The probability of error is also invariant to transformations in which 
each of the signals is right-multiplied by a respective MXM unitary matrix. There may be a distinct such matrix for each 
signal in the constellation. (Right-multiplication is sometimes referred to as "post-multiplication.") 
[0054] Below, we will describe a novel method for constructing a signal constellation, variations of such a constellation 
are readily obtained by these transformations involving left- or right-multiplication by unitary matrices. Such variations 
are to be regarded as equivalent to the constellations obtained by the method described below, and thus to lie fully 
within the scope of our invention. 

[0055] Turning to FIG. 2, our method begins at block 25 with the step of providing an initial signal matrix 0> v Any 
TxM matrix is acceptable for this purpose, provided it is otherwise a valid signal matrix; that is, provided its columns 
are all mutually orthonormal. (An alternative way to state this condition is that = I , where I is the MxM identity 





<*=max|K*rl • 
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matrix.) 

[0056] As noted above, a set of L signal matrices is readily constructed from the initial matrix and a TXT unitary 
matrix, which is referred to as the generator matrix. The initial matrix is left-multiplied by the generator matrix p times 
p ~ 1> 2 L " 1 ' ,0 9 enera,e '"rther signal matrices of the constellation. A set of size 4, larger than L may first be 
generated and then a constellation of size L specified by selecting a proper subset of the generated signal matrices 
This selection is advantageously made in such a way that signals that tend to reduce the error probability are favored 
and signals that tend to increase the error probability are disfavored. 

[0057] In the discussion that follows, it is assumed, for simplicity of presentation, that L is equal to the full number 

I n 9 en °f ed rf s, 9 nal ma,rices Extensions of this discussion to cases in which L < J will be readily apparent to those 
SKHI6CJ in in© art. 

[0058] One class of generator matrices that we have found particularly advantageous is the class of matrices that 
are L th roots of the TxT identity matrix l T . For convenience (and without loss of generality), it may be assumed that 
each of these matrices is a diagonal matrix; that is. the T diagonal elements are finite, and all other elements are zero 
Each of the d.agonal elements is a complex number whose absolute value is 1. (As used herein, the term 'complex 
number means a number that is pure real, pure imaginary, or a sum of real and imaginary parts.) Each of these 
elements is expressed in the form 



where e is the root of the natural logarithms, r is the row (or column) index of the pertinent matrix element (1=12 
T). and u, is ; an integer on the closed interval from 0 to L-1. (In fact, a more general class of generator matrices is 
readily constructed by letting the factors u,take on continuous, rather than integral, values ) 

[0059] It will be appreciated that when such a generator matrix is repeatedly applied as a left-multiplier of the initial 
thTin'tiarmatrix S '' 9nal maM ° eS ° an be 9 enerated However, the Lth application will simply reproduce 

[0060] Those skilled in the art will also appreciate that given any pair of distinct signal matrices, the corresponding 
correlation matrix will be invariant to left-multiplication of these signal matrices by the generator matrix This will be 
™£« anV " Umber ° f re P eated "eft-multiplications, provided the number is the same for both signal matrices 
hSSi AS 3 consec ' uence - th e constellation of L generated matrices has a circulant correlation structure 
[0062] As a further consequence, given the pair <t>„ the generator matrix can, in particular, be applied just enough 
times to transform matrix 0, into matrix <&, . Therefore, the set of correlation matrices *,t$, fc2 3 L exhausts 
all possibilities (except the self-correlation, which is not of interest here). 

[0063] If a generator matrix 0 is an Lth root of l T . ft is completely specified by the vector u = (u, u T ) This 

vector, in turn, completely specifies the corresponding signal constellation of size L generated from a given initial signal 
vector. In general, there are very many distinct possible valuations for the components of this vector It is desirable to 
select a valuation that tends to minimize the error probability. 

[0064] One approach is to perform an exhaustive search, over all possible valuations of u, for that u-vector that yields 
the smallest value of, e.g., 5 (for a given initial signal matrix). However, in many practical situations, Land T will be so 
nn e e c II 6 arS m ° re POSSib,e -y- vectors tnan can be searched in an economically reasonable amount of time 
on Hn L o 6 « ' ° ne US6fUl alternate approach is to generate a random selection of u-vectors, as indicated at block 
30 of FIG. 2. As ind.cated in cycle 35 of the figure, and more specifically by block 40, a constellation of L signal matrices 
•s generated, in turn, from the generator matrix represented by each of the u-vectors. As indicated in sub-cycle 45 
and more specifically by block 50, an exhaustive set of L-1 correlation matrices is formed for each of these constella- 
tions. As md.cated by block 55. a measure of error probability is computed for each of these correlation matrices One 
useful measure is the C.U. B. Another useful measure is the norm Btft^j), which we denote in the figure by the symbol 8, 

(Thus. 0 = max 6, .) 



t<ISL 



[0066] After sub-cycle 45, but still within cycle 35, the greatest of all of these error-probability measures is identified 
(for the current constellation), as indicated by block 60. After this maximum error-probability measure has been deter- 
mined for each constellation (and thus, after cycle 35 has been completed for all u-vectors). the optimal constellation 
is chosen. As md.cated by block 65. the chosen constellation is that constellation whose maximum measure of error 
probability is smallest. 

[0067] A new constellation can be generated by either changing the initial signal matrix 0>, or changing the u-vector 
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(which is equivalent to changing the generator matrix). The procedure of FIG. 2 is a one-parameter search, in which 
only the u-vectors were varied. An alternative is the two-parameter search procedure of FIG. 3 : in which both u-vectors 
and the initial matrix are varied. (However, we have found that in general, performance is more sensitive to the choice 
of u than it is to the choice of O v ) 

5 [0068] Specific techniques for carrying out two-parameter searches are well known in the art, and need not be de- 
scribed here in detail. Accordingly, FIG. 3 illustrates an exemplary two-parameter search at a high level. With reference 
to the figure, a set of initial matrices is generated at block 70, and a set of random u-vectors is generated, as before, 
at block 75. At block 80, a constellation of signal matrices is generated for each possible pairing of an initial matrix with 
a u-vector. At block 85, the maximum error-probability measure is determined for each of these constellations. At block 

10 90, an optimal constellation is chosen. As before, the optimal constellation is that constellation whose maximum error- 
probability measure is the smallest. 

[0069] It should be noted that the precise blocks illustrated in FIGS. 2 and 3 are intended as a pedagogical aide. 
Those skilled in the art will appreciate that procedures involving computational shortcuts and other departures from 
the illustrated flowcharts can lead to the same overall results. AH such procedures should be deemed to be equivalents 
is of the illustrated procedures and to lie within the scope and spirit of the invention. 

[0070] It should be noted that not all possible initial signal matrices O n are equally good. For example, one undesirable 
matrix is the matrix whose //-th element (i is the row index and /is the column index) is 1 if / = /, and 0 otherwise. This 
matrix will lead to correlation matrices all of whose singular values are 1. 

[0071] The matrix described above has all zeroes for the entries in the last T-M rows. In fact, we believe that it is 
20 also disadvantageous to use for the initial signal matrix any matrix having T - M rows of all zeroes. 

[0072] We now describe one useful method that we have found for constructing an initial signal matrix. This method 
gives a relatively high probability of generating a good initial signal matrix. 

[0073] A Discrete Fourier Transform (DFT) matrix is a square matrix, having T orthonormal columns, in which the 
rs-th element (r is the row index and s is the column index) is given by 

25 

^(M)(S-1) 1 
e X~p. 

30 According to our method, the initial signal matrix is constructed by selecting M of the T columns of a DFT matrix. Such 
a selection may be made randomly, for example. 

[0074] According to an alternate method, the initial signal matrix is constructed as a random-valued orthonormal 
matrix. For example, procedures are readily available for generating a TxM matrix X whose elements are T-M inde- 
pendent complex variables, Gaussian-distributed with 0 mean and unit variance. A well-known procedure referred to 
35 as "QR factorization" is performed on this matrix. QR factorization may be represented by the equation X = QR, where 
Q is a TxM matrix whose columns are orthonormal, and R is an upper or lower-triangular square matrix. The initial 
signal matrix is set equal to Q. 

[0075] If, for example, a signal constellation or an initial signal matrix is newly chosen or newly generated at the 
transmitting location, information sufficient to specify such constellation or matrix is then advantageously transmitted 

40 to the receiving location for use in the receiver. Conversely, the signal constellation or initial signal matrix may be 
chosen or generated at the receiving location, and transmitted to the sending location for use in the transmitter. 
[0076] According to one extension of the ideas described above, multiple generator matrices may be used in a 
composite fashion to generate a signal constellation from an intitial signal matrix. That is, let the matrix G-, be an Ljth 
root of l T , and let the matrix 0 2 be an l^th root of l T . Then a constellation is readily generated by applying, to an initial 

45 signal matrix, compound rotations of the form e^efe" 1 , /, = 1,..., L, ; ^ =1, .... . The resulting constellation will 
have a circulant correlation structure described by the^ollowing rule: 

Given indices c^, ^ f v p 2 , y v y 2 , 5,, 8^ if (o^ - fc) mod L n = ft, - 5^ mod L,, and (o^ - p 2 ) mod l_2 = (? 2 - ) 
mod Lg, then Ot aia2 O pi p 2 = 0t Tm o 6l62 . 

[0077] Those skilled in the art will appreciate that this technique can be extended to a composite of an arbitrary 
so number K of generator matrices of the relevant kind, so that each generator matrix is an L' k th root of an identity matrix, 
k=1,..., K. Thus, eg., L, x «. x Impossible USTM matrices may be generated. 

[0078] Thus far, we have described procedures that are particularly appropriate when the characteristics of the phys- 
ical transmission channel between the transmitting and receiving antennas (i.e., the elements of the matrix H) are 
unknown, or are assumed to be unknown. With certain modifications that are to be described, these procedures are 
55 also useful when the channel characteristics are known. 

[0079] When the elements of the matrix H are known, the maximum likelihood decoding is given by: 
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= arg nun tr { (X - V^cD, H) (X - J%<P, H)' } . 
s The pertinent Chernoff upper bound is given by: 

J A 1 



c.u.B.=^n F 

10 m =\ 1 
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In the preceding expression, the quantities D m are the singular values of difference matrices of the form <t>, - 1* As 
cor?eTa«on% S trS e re * 6Xamine diff6renCe ma,riC6S (e ' 9 - ' ** '= 2 L >■ because there remains a circulant 

E am* tl dent ,ha L n ° W - the P U B iS minimized «hen all of the singular values are driven as high as possible 
[0081] All of the procedures described above for generating u-vectors and initial signal matrices also apply in the 
present, known-channel case. However, the optimization procedures for finding the best possible signal constellations 
are now directed to maximizing the singular values of difference matrices, and not to minimizing the singular values 
of correlation matrices. 

[0082] In practice, the signal constellation may be stored at the site where signals are generated for transmission 
or at a location accessible to it. Alternatively, the signal matrices may be generated, as needed, at such a site or 
location. Generally, it will be advantageous to generate the signal matrices as needed, because much less memory 
will be taken up in that case. That is, instead of storing a large number of signal matrices, it will suffice to store information 
hat defines the initial signal matrix and the generator matrix or matrices. (If the L signal matrices in the signal constel- 
lation are fewer than the J generated matrices, it will also generally be necessary to store a rule indicating how the L 
signal matrices are to be chosen.) 

[0083] As noted the signal matrices are also used, at the receiving end, for recovering the transmitted signal from 
the received signal. For this purpose, again, the signal constellation may be stored where it will be accessible or it 
may be generated as needed. Generally, it will be advantageous to generate the signal matrices as needed for the 
same reasons discussed above. ' 
[0084] A USTM signal constellation of certain specific kinds we have described here may be characterized in various 
ways, e.g., as a subset of a set of signal matrices having a circulant correlation structure, or as a set of signal matrices 
Z k 10 w •u t 1 er Via re P eated multiplications of a common initial signal matrix by one or more unitary matrices 
We have described a particular mathematical formalism, specifically a matrix formalism, for generating these signal 
constellations. It should be noted that other mathematical formalisms, such as formalisms related to the theory of seff- 
correctmg codes, may be employed to achieve an equivalent end result. Thus, the scope of our invention is not defined 
by the specific mathematical formalism that is invoked for generating the signal constellation, but instead is meant to 
encompass all such formalisms.provided equivalent end results are achieved. 

Example 

[0085] We generated a constellation of £ = 257 signal matrices having T = 8 and M = 3. by methods substantially 
as described above. We assumed that the elements of the matrix H (i.e.. the propagation coefficients) are unknown 
The u-vector from which this constellation was generated was: 

(220 191 6 87 219 236 173 170). The initial signal matrix was: 



55 
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The error-probability measure was evaluated as 5 = 0.74355150. 

[0086] We carried out simulation experiments using the first L=256 signals of this constellation, in the unknown- 
channel case. FIG. 4 shows the block error rate (i.e., the probability of falsely identifying a received signal matrix), and 
FIG. 5 shows the bit error rate that we obtained, both at a data rate of one bit per channel use. It is evident that the bit 
error rate drops steeply as the signal-to-noise ratio p is increased. This is consistent with our theoretical prediction that 
USTM signals are asymptotically optimal for large signal-to-noise ratios. 



Claims 

1. A method of wireless communication that comprises generating a plurality of baseband signals, modulating the 
signals on a carrier, thereby to form a sequence of carrier-level signals, and transmitting the carrier-level signal 
sequence from an array of one or more antennas, wherein: 

each baseband signal comprises one or more sequences, in time, of complex amplitudes, each said sequence 
to be transmitted from a respective antenna of the array; 

each baseband signal is representable as a transmission matrix in which each column represents a respective 
antenna and each row represents a respective time segment; 

each transmission matrix is proportional to a respective matrix having orthonormal columns, said matrices 
having orthononnal columns to be referred to as USTM matrices; and 

the generating step comprises generating a plurality of USTM matrices, said generating carried out such that 
each resulting USTM matrix is either a common initial matrix or a product of one or more multiplications of the 
common initial matrix by one or more unitary matrices. 

2. The method of claim 1 , wherein each baseband signal is proportional to one of the generated USTM matrices. 

3. The method of claim 1 , wherein each baseband signal is proportional to a matrix derived from one of the generated 
USTM matrices by multiplication by a unitary matrix. 

4. The method of claim 1 , wherein the baseband signals are made from some, but not all, of L possible USTM matrices 
that may be generated by the generator matrix. 

5. The method of claim 1, wherein the baseband signals are made from some, but not all, of L, x ... x L K possible 
USTM matrices that may be generated by the generator matrices. 

6. The method of claim 1 , further comprising transmitting a description of selected common initial matrix to a receiver. 

7. A method of wireless communication that comprises generating a plurality of baseband signals, modulating the 
signals on a carrier, thereby to form a sequence of carrier-level signals, and transmitting the carrier-level signal 
sequence from an array of one or more antennas, wherein: 
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to^e urnfmS2?r^? mPriS ^ ^ ° F ™** se ^ uences * in { ™ e > of com P J ** amplitudes, each said sequence 
to De transmitted from a respective antenna of the array 

ISIESS^ 38 3 ,ranSmission ma,rix in *** -lumn represents a respective 

anienna ana each row represents a respective time segment 

^rnn r oT iSSi ° n , ma ! riX iS proportional to a respective mat'rix having orthonormal columns, said matrices 
having orthonormal columns to be referred to as USTM matrices 

lesuSSus al P bseTof ri r, S ? h'^ 0 * "f °' US ™ ma,riCeS ' Sa ' d Carried oul «•* *- the 

resulting set is a subset of a set having a circulant correlation structure and 

each transmission matrix is proportional to a USTM matrix from the subset, or to an equivalent of such a matrix 
derived from such a matrix by multiplication by a unitary matrix. 

ln^ of ^ ireles = com v munica «on that comprises generating a plurality of baseband signals, modulating the 

8 S6qUenCe ° f CarrieMeVe ' Si9na,S ' and tranSmittin9 ,he ca-ieMeve sfgna! 
sequence irom an array of one or more antennas, wherein: 

tobe SnsSedC^r 6 , 8 ^ 7 SeqUenC6S ' in ^ ° f COmp,6X am P'«^es, each said sequence 
io De transmitted from a respective antenna of the array 

ZennTZnT '* fepreSen,ab,e as a transmission matrix in which each column represents a respective 
antenna and each row represents a respective time segment ^ 

each transmission matrix is proportional to a respective matrix having orthonormal columns said matrices 
having orthonormal columns to be referred to as USTM matrices- matrices 
the generating step comprises selecting USTM matrices from a set of USTM matrices in which each USTM 

In 6 !^ ° f WifeleS !,: communication that comprises generating a plurality of baseband signals, modulating the 

IZ^e^ZaZTJl f ° rm 3 S T enee °' CarrieNeVe ' Si9na ' S ' 3nd ,ranSmi,,in9 ,ha ca'rrieHetrsfgna 
sequence irom an array of one or more antennas, wherein: 

tobe ^m^nf 31 C ° mprises ° ne or more sequences, in time, of complex amplitudes, each said sequence 
to be transmitted from a respective antenna of the array 

each baseband signal is representable as a transmission matrix in which each column represents a respective 
antenna and each row represents a respective time segment- represents a respective 

hav^ r oZ™, ma | riX iS P t T T ni0na * ,0 3 r6SpeCtiVe ma, ' rix havin 9 °nhonormal columns, said matrices 
having orthonormal columns to be referred to as USTM matrices- 

%£Z£S2S> ZZSZST 0 us ™ ma,rices ,rom a subset of a set of us ™ ^ a 

each transmission matrix is proportional to one of the selected USTM matrices, or to an equivalent of such a 
matrix derived from such a matrix by multiplication by a unitary matrix. equ.va.ent of such a 

10. A method of wireless communication that comprises receiving a carrier-level signal sequence from an array of one 

^sssssFxsr to a transmission ' and de ™ 9 said ™ ^ ^ 

^ele^ 

each baseband signal is representable as a matrix in which each column represents a respective antenna and 
each row represents a respective time segment respective anienna ana 

£«oT et f h ° d f T ef COmprises ascribin 9 to eac * ^ceived matrix a signal matrix belonging to a known constel- 
at.on of signal matrices, to be referred to as USTM matrices, having orthonormal coLns- 

^Tr^ZT n ma , triX, K hS aSCribinQ StSP comprises comparing the received matrix to each USTM matrix 
of the constellation, thereby to evaluate a likelihood score, and selecting for ascription that USTM matrix that 
yields the greatest likelihood score; and 

laXt^i^ 9 r ratin9 3 PlUraHty ° f US ™ ma,riCeS ' Said 9eneratin 9 camed ° ut ^ch that 

common^ Z ~ PJ T " * C ° mm ° n M ma,riX ° r a product * one or more multiplications of the 
common initial matrix by one or more unitary matrices. 
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11. The method of claim 7, wherein each USTM matrix in the constellation is one ol the generated USTM matrices. 

12. The methodpf claim 7, wherein each USTM matrix in the constellation is derived from one of the generated USTM 
matrices by multiplication by a unitary matrix. 

13. The method of claim 1 or 7 wherein only one unitary matrix is used to generate USTM matrices by repeated 
multiplications, said unitary matrix to be referred to as a generator matrix. 

14. The method of claim 13, wherein the generator matrix is an/ th root of an identity matrix, / a positive integer. 

15. The method of claim 11, wherein the constellation contains some, but not all, ol^f possible USTM matrices that 
may be generated by the generator matrix. 

16. The method of claim 1 or 7 wherein two or more unitary matrices are used to generate USTM matrices by repeated 
15 multiplications, said unitary matrices to be referred to as generator matrices. 

17. The method of claim 13 or 16 wherein: there are K generator matrices, K a positive integer greater than 1; each 
generator matrix is a-4b'th root of an identity matrix, k= 1...., K; and each-4 is distinct. 

20 18. The method of claim 1 4, wherein the constellation contains some, but not all, o\A X-xA possible USTM matrices 
that may be generated by the generator matrices. 

19. The method of claim 1 or 7 further comprising selecting the common initial matrix. 

25 20. The method of claim 1 9, wherein the common matrix selection is carried out by selecting columns of a DFT matrix. 

21 . The method of claim 1 9, wherein the common matrix selection is carried out by generating a matrix having random- 
valued complex coefficients and performing a QR factorization on said matrix. 

30 22. The method of claim 19, further comprising transmitting a description of the selected common initial matrix to a 
station whose signals are to be received. 

23. A method of wireless communication that comprises receiving a carrier-level signal sequence from an array of one 
or more receiving antennas in response to a transmission, and demodulating said sequence, thereby to form a 

35 plurality of baseband signals, wherein: 

each baseband signal comprises one or more sequences, in time, of complex amplitudes, each said sequence 
received from a respective antenna of the array; 

each baseband signal is representable as a matrix in which each column represents a respective antenna and 

40 each row represents a respective time segment; 

the method further comprises ascribing to each received matrix a signal matrix belonging to a known constel- 
lation of signal matrices, to be referred to as USTM matrices, having orthonormal columns; 
for each received matrix, the ascribing step comprises comparing the received matrix to each USTM matrix 
of the constellation, thereby to evaluate a likelihood score, and selecting for ascription that USTM matrix that 

45 yields the greatest likelihood score; and 

the method further comprises generating a set of USTM matrices, said generating carried out such that the 
resulting set has a circulant correlation structure. 

24. A method of wireless communication that comprises receiving a carrier-level signal sequence from an array of one 
so or more receiving antennas in response to a transmission, and demodulating said sequence, thereby to form a 

plurality of baseband signals, wherein: 

each baseband signal comprises one or more sequences, in time, of complex amplitudes, each said sequence 
received from a respective antenna of the array; 
55 each baseband signal is representable as a matrix in which each column represents a respective antenna and 

each row represents a respective time segment; 

the method further comprises ascribing to each received matrix a signal matrix belonging to a known constel- 
lation of signal matrices, to be referred to as USTM matrices, having orthonormal columns; 
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for each received matrix, the ascribing step comprises comparing the received matrix to each USTM matrix 
of the constellation, thereby to evaluate a likelihood score, and selecting for ascription that USTM matrix that 
yields the greatest likelihood score; and 

each ijSTM matrix in the constellation is either a common initial matrix, or a product of one or more multipli- 
cations of the common initial matrix by one or more unitary matrices. 

25. A method of wireless communication that comprises receiving a carrier-level signal sequence from an array of one 
or more receiving antennas in response to a transmission, and demodulating said sequence, thereby to form a 
plurality of baseband signals, wherein: 

each baseband signal comprises one or more sequences, in time, of complex amplitudes, each said sequence 
received from a respective antenna of the array; 

each baseband signal is representable as a matrix in which each column represents a respective antenna and 
each row represents a respective time segment; 

the method further comprises ascribing to each received matrix a signal matrix belonging to a known constel- 
lation of signal matrices, to be referred to as USTM matrices, having orthonormal columns; 
for each received matrix, the ascribing step comprises comparing the received matrix to each USTM matrix 
of the constellation, thereby to evaluate a likelihood score, and selecting for ascription that USTM matrix that 
yields the greatest likelihood score; and 

the constellation is a subset of a set of USTM matrices having a circulant correlation structure. 
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